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Abstract 

We apply the functional renormalisation group to few-nucleon systems. Our starting point is 
a local effective action that includes three- and four-nucleon interactions, expressed in terms of 
nucleon and two-nucleon boson fields. The evolution of the coupling constants in this action is 
described by a renormalisation group flow. We derive these flow equations both in the limit of 
exact Wigner SU(4) symmetry and in the realistic case of broken symmetry. In the symmetric 
limit we find that the renormalisation flow equations decouple, and can be combined into two 
sets, one of which matches the known results for bosons, and the other result matches the one for 
fermions with spin degrees only. The equations show universal features in the unitary limit, which 
is obtained when the two-body scattering length tends to infinity. We calculate the spin-quartet 
neutron-deuteron scattering length and the deuteron-deuteron scattering lengths in the spin-singlet 
and quintet channels. 

PACS numbers: 21.45.-v Few-body systems, 24.10.Eq Coupled-channel and distorted-wave models, 25.45.-z 
2H- induced reactions, 11. 10. Hi Renormalisation group evolution of parameters 

Keywords: Deuteron-deuteron scattering, nucleon-deuteron scattering, functional renormalisation group, 
SU(4) symmetry, Wigner supcrmultiplet 
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I. INTRODUCTION 



Few-body physics provides a solid starting point for analysing various nonperturbative 
approaches and methods which can then be used for studying more complicated many-body 
problems. The study of few-body systems also allow us to determine the input parameters 
needed for those complex problems, without the need to rely on further approximations. 
One promising method is the Functional Renormalisation Group (FRG) l|, |2j], as it provides 
a framework that can simultaneously describe both few- and many-body systems. Some 
reviews of the method and its applications can be found in Refs. js-5]. 

The approach is based on a running effective action (REA) which is a generalisation of 
the standard quantum mechanical effective action - the generating functional of the one- 
particle irreducible Green functions. The REA includes the effects of all the fluctuations with 
momenta in the region q 2 > k 2 where k is a running scale. This is achieved by introducing 
a regulator that, either sharply or smoothly, suppresses the contribution of modes with 
q 2 < k 2 . The FRG describes the evolution of this action as the cut-off scale k is lowered. As 
k approaches zero, all fluctuations are included and thus the full effective action is recovered. 

The FRG method has already been applied to few-body problems in a number of papers 
-|l0|. In Ref. jfj], it was exploited to derive the Skornyakov-Ter-Martirosyan equation 
for three-body systems [ll]. The Efimov effect 
addressed in Ref 



12] in bosonic and fermionic systems was 
f. An important extension of the approach was developed in Ref. js] 



to treat the four-body problem in the presence of the three-body Efimov effect, and the 



process of dimer-dimer scattering was studied in detail in Ref. 



10] . This shows one of the 



key strengths of the FRG method, which is its universality. With relatively minor changes 



it can be adopted to study a large variety of problems in particle, nuc 



13 



ear and condensed 



4l5|. 



matter physics. Some representative examples can be found in Refs. 

Up to now, most applications of the FRG have been to particle physics and condensed 
matter systems. In the present work we apply it to systems of up to four nucleons, laying 
the groundwork for extensions to larger numbers of nucleons and to nuclear matter. One 
important aspect of our implementation of the FRG is the introduction of bosonic fields to 
describe interacting pairs of nucleons. This will be very useful for future work since it will 
allow us to link the scattering of nucleons in vacuum to pairing in systems of many nucleons, 



as in Ref. 



13] . This study of few-nucleon systems allows us to develop some of the tools and 



determine the parameters that will be needed in such work. In this context, it is helpful to 



make use of Wigner's SU(4) "supermultip 



approximate symmetry for light nuclei 



16]. 



et" symmetry, which is known to be a very good 



The running effective action, denoted by T, satisfies an exact flow equation, a functional 
differential equation of the form 

d k T = -- Str [(d k R) (T^ - Ry 1 } . (1) 

where r^ 2 ^ denotes the second functional derivative taken with respect to the fields enter- 
ing the action, and R is a matrix of cut-off functions. These functions act as regulators, 
suppressing the contributions of fluctuations with momenta below the running scale k, and 
thus driving the evolution of the system as k is lowered. The structure of the evolution 
equation is rather straightforward, the main complexity is in the operation Str, which is the 
supertrace, which is taken over both energy-momentum variables and internal indices. It is 
needed since we consider a mixed system of fermions (nucleons) and bosonic dimers in this 
work. 

In practical applications, the REA is usually truncated to a finite number of local terms. 
The choice of these is guided by the relevant physics in the system we wish to describe. This 
truncation reduces the full functional differential equation ([I]) to a set of coupled ODE's for 
the running coupling constants and renormalisation factors multiplying those terms in the 
action. These equations are straightforward to solve using standard methods. 

The use of local interactions means that the FRG method has close links to approaches 



based on effective field theories 13, llSj • In some cases, for example the action used to study 
the Efimov effect in three-boson systems, the terms included are just the leading terms of the 
corresponding effective field theory. In other cases, such as the applications to dense matter, 
strict power counting arguments cannot be used and we must rely on the need to describe 
emergent aspects of the physics that are known to be important, such as superfluidity. 

There is natural scaling of the evolution equations for large k, where the results coincide 
with the behaviour of the equations in the unitary limit, when the scattering length is 
infinite. In that limit we get simple universal equations. Thus all the evolution equations 
must collapse to these universal ones as we increase the cut-off scale k. The universal 
equations can then be used to obtain boundary conditions for the evolution equations away 
from the unitary limit. 



In the main body of the paper we derive a general set of equations for the evolution of the 
REA, without assuming SU(4) symmetry. The special case of exact symmetry is analysed 
in detail in the Appendix, where we show that there are two decoupled sets of equations. 

One of these sets describes the channels where the nucleons are in states of mixed spatial 
symmetry and it has the same form as that for a system of fermions with spin degrees 
of freedom only, as derived in Ref. 10j. The other set describes the spatially symmetric 
channels and, after a suitable redefinition of couplingconstants, it matches the evolution 
equations for a system of interacting bosons in Ref. 7|. Despite the similar forms of the 
two sets of equations, the differences in their numerical coefficients have profound physical 
consequences. In particular, the equations for spatially symmetric systems leads to the 



Efimov effect 12[ and limit-cycle behaviour of the 3-body couplings 19[. In contrast these 



effects do not occur for few-fermion systems without spatial symmetry. 

In the next section we set out the form of the REA we use in our applications of the FRG 
to few-nucleon systems, and in Sec. Ill we obtain the evolution equations for the running 
parameters in that action. Then, in Section IV, we apply these to analyse scattering in three- 
and four-nucleon systems for realistic NN scattering lengths. We compare our results to 
those in the SU(4) limit, details of which can be found in Appendix A. Finally, some technical 
details of the link between scattering lengths and REA are discussed in Appendix B. 



II. RUNNING EFFECTIVE ACTION 

Since our ultimate goal is to study pairing in dense nuclear matter, we must isolate the 
relevant degrees of freedom. As is well known, the nuclear force is just strong enough to 
generate a bound state, the deuteron, in one of the S"-wave channels. Its isospin-1 analogue 
is only just unbound and appears as a virtual state very close to threshold. We thus find 
it useful to introduce boson fields to describe the lowest two-body states in both these 
channels. If one imagines starting from an effective action with a local two-body interaction 
between the fermions, the boson fields can be introduced using a Hubbard-Stratonovitch 
transformation, which is exact in this case. We use the notation t for the (spin-triplet) 
vector-isoscalar boson field, corresponding to the deuteron, and s for the (spin- singlet) 
scalar-isovector one. These we refer to collectively as "dimers", and label by a capital D. 
We reserve the lower-case label d for a deuteron, and N for a generic nucleon. Finally, when 
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specifying quantum numbers of particular channels, we use the 5*0(4) notation (S, I), where 
S and / denote the resultant spin and isospin quantum numbers, respectively. In the present 
we work, we assume that isospin is a good symmetry. 

The key approximation we make is to truncate the full REA to one with only a finite 
number of parameters, so that the functional differential equation can be replaced by a set 
of ordinary differential equations. In the present case, we introduce a set of local contact 
interactions in the channels of interest, as well as kinetic terms for the dimer fields. The 
form for the REA we use is 

P 2 



+ j d 4 pt\(p ,p) (z^tpo - Z m,t-^ - u 1)t + ie \ U(p ,p) 
+ j d A ps\(p 0) p) (z^sPq - z m, s j^ - ui t , + iej s a (p ,p) 

+r 2 + r 3 + r 4 . (2) 

Since the action at the start of the evolution is obtained by bosonising a purely fermionic 
one, the boson fields should be non-propagating auxiliary fields for k — > oo. Their wave 
function and kinetic-mass renormalisation factors {Z^ and Z m ^) should tend to zero in this 
limit. When we run the action as we lower the cutoff scale k, these fields become dynamical, 
and the renormalisation factors grow. 

The bosonisation introduces self-energies for the boson fields (/ui,j) and couplings of the 
dimers to pairs of nucleons in the (1,0) and (0, 1) two-body channels. These couplings are 
described by the term, 

r 2 = j 5 (4) (pi +p 2 - P3 )d' l p 1 d'p 2 d 4 p 3 V2 (3) 

V 2 = -^fo9t * f (P30,P 3 ) • [^(Pio,Pi)^ C (P2o,P 2 )] (1 ' 0) + h -c- 

1 r + , . r , . , ,n, ,i(0,l) 



s\P30,P 3 ) ■ [^(Pio,Pi)^ C (P2o,P 2 )] ' +h.c. , (4) 



where 



' l l ) mt 2 m S2 T rn Bl m S2 < J mt 1 mt 2 ' l l J ™-t 1 m Bl - (5) 

Since is a (1,0) tensor and a (0,1) one, the scalar products in Eq. fl3j) should be 
understood as being taken in either spin or isospin space, as appropriate. The values of the 
coupling constants depend on a choice of scale in the dimer fields when they are introduced 



by the Hubbard-Stratonovitch transformation. Thus we must find that the values of the gi 
do not appear separately in physical quantities, but only in combinations such as g\ju\^. 
An illustration of this is provided by the discussion of the evolution of the in the next 
section. The unobservable coupling constants g^ do not run in vacuum and, for simplicity, 
we choose their values to be equal, 

9t,s = 9- (6) 

To describe the three-nucleon channels, we introduce a set of local dimer-nucleon inter- 
actions, described by the term 

r 3 = - J 5 (4) (pi +P3-P2- P4)d 4 pid 4 p2d 4 p 3 d 4 p 4 V 3 , DN (7) 

i,j=t,s 

+AL 3Al/2) [*W^WPi)] (3/2 ' 1/2) • [^ 2 o,P 2 )^4o,P 4 )] (3/2 ' 1/2) 
+AE /2 ' 3/2) [ S t (P30,P 3 )^(Pio,P 1 )] (1A3/2) • [^(p 2 o,P 2 )s(p 4 o,P 4 )] (1/2 ' 3/2) , (8) 

where A^ 2,1 ^ 2 '' = X^/ 2 ' 1 ^ 2 ^. These have been expressed in terms of interactions in the 
doublet-doublet channel, with spin-isospin quantum numbers (1/2,1/2), and the quartet- 
doublet channels, with quantum numbers (3/2,1/2) or (1/2,3/2). The (1/2,1/2) channel 
has the quantum numbers of the ground states of 3 H and 3 He. 

Finally, we introduce one class of four-nucleon interactions, represented by local two-body 
dimer-dimer interactions and described by the term 

T 2 ,dd = - J <5 (4) (pi +P3-P2- P4>d 4 p 1 d 4 p 2 d' l p3d 4 p i V2,DD, (9) 

V 2 ,DD = -^2u 2 , i (i l! ip W ,P 1 ) ■ «(p20,P 2 )* t (P30,P3) ' *(P40,P 4 ) 

i=t,s ^ 

-g* + OlO,Pl) •* t (P30,P 3 )*(P20,P 2 ) -i(P40iP4, 



+U 2 ,ts\t' t (p 1 o,Pi) ■ t{p W ,P 2 ) S f b30,P 3 ) -S(P40,P 4 ) 

+ E U2,iji\piO,Pl) ■ i\p30,P 3 ) j(P20,P 2 ) ■ J'(P40,P 4 ) ; (10) 

ij=t,s 

where u 2 ^ s = u 2 , s t- The u 2 ,i terms act in the (2, 0) and (0, 2) two-dimer (four-nucleon) 
channels, the u 2 j s terms act in the (1, 1) channel, and the u 2 terms act in the (0, 0) channel. 
Only the last of these has the quantum numbers of the ground state of the a particle. 
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The behaviour for large k is particularly simple. In this region, the regulator provides 
the only scale and every coupling varies as a power of k. The effects of SU(4)-breaking 
also become negligible here and so the evolution equations reduce to the simpler, SU(4)- 
symmetric forms discussed in Appendix A. This allows us to determine the initial conditions 
on the evolution here and these then provide the physical input into our approach. 

Most of the parameters are "irrelevant" in the technical sense that the low-energy results 
from our action in the physical limit (k — > 0) do not depend strongly on the precise values 
of the initial conditions imposed at very large k. However, there are three exceptions and 

re self-energies, Un, 



2l|. Their values at 



these provide the physical input into our approach. Two of them are t 
which are "relevant" parameters in RG language, as discussed in Ref. 
k = determine the strength of the interaction between the fermions and can be related 
either to the NN scattering lengths, a t and a s , or, in the channel with the bound state, to 
the deuteron binding energy, 

fc = -T^- (11) 



Ma 



The third parameter is associated with the spatially-symmetric channel of the three- 
nucleon system. This is the channel that displays the Efimov effect 12] and as a result 
the evolution of its coupling constant for it has a limit-cycle behaviour at large- /c. Just 
as in the analogous EFT treatment 19j, one piece of three-body data is needed to fix the 
starting point on this cycle. Here we choose the nucleon-deuteron scattering length in the 
spin-doublet channel, and explore how four-nucleon observables are related to it. 



III. EVOLUTION EQUATIONS 

To derive the evolution equations for the coupling constants in our action, we substi- 
tute the parametrisation given in the previous section into the right-hand side of Eq. (JTj). 
The inverse of the second derivative of T forms a scale-dependent propagator, and the driv- 
ing terms for the coupling constants all have structure of one-loop integrals, with differing 
numbers of external fields. In general these include expressions with powers of fields or 
derivatives beyond those contained in our truncated action. We therefore need to expand 
the loop integrals in powers of the fields and non-localities (i.e., in energies and momenta) 
and pick out those terms that match the structures listed in the previous section. The 
scale-dependent coefficient in front of each structure then gives us the driving term in the 
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differential equation for the corresponding coupling constant. 

As well as the choice of truncation, we also need to make a choice of the point around 
which we expand our energies. For the dimers, we choose the energy of the lowest two- 
body bound state, In the four-nucleon sector, for example, this means that we expand 
around the deuteron-deuteron threshold. For the nucleons we expand around one half of 
this binding energy. Thus in the nuclear-dimer channels of the the three-nucleon system we 
expand about an energy Sd/2 below the the scattering threshold. This procedure can be 
justified rigorously in the limit of exact SU(4) symmetry, as discussed further below, but 
it is only approximate in the general case. It could thus lead to a slower convergence than 
expected of the expansion. 

Finally, we need to choose the forms for the cut-off functions in the regulators for the 
nucleons and dimers. Here we take the form suggested by Litim 



suited for partially analytic calculations. It is also optimised for actions truncated to purely 
local (energy- independent) interactions, as also discussed in more detail by Pawlowski 
The regulators are thus given by 



20| for both, as it is well- 
irely 



RN(q;k) = ^fe(k-q), (12) 
R D (q;k) = Z^k)-^e{k-q). (13) 

Note that we have assumed that the two Z^s are equal; this can be shown to hold in 
the expansion scheme we use here, see Eq. (|23|) below. The inclusion of the wave-function 
renormalisation factor in the definition of Rd has the advantage of allowing us to scale out 
the parameters g, M and at, leaving much simpler, dimensionless expressions. For example, 
a generic three-body coupling A has a natural scaling of the form 

AW = ^ A( *>- < 14) 

Analogously a generic four-body coupling u 2 can be scaled as 

U 2 {k) = g- 4 M- 3 at 3 u 2 (k). (15) 

A. One-boson terms 

The evolution equations for un, Z^i and Z m ^ are relatively straightforward in vacuum. 



Their forms are the same as in the simpler case studied in Ref. 



101 ]; For example, both the 



U\i satisfy the differential equation 



OkU lt i{k) = —-—3 / drq- -j, (16) 



2 (2nY 



2M 



R N (q; k)-S d /2 



where Ed appears in the denominator as a result of our choice of expansion point and we 
have already performed the integral over the virtual energy q°. This equation is an exact 
differential and can be integrated directly. Its solution grows linearly with k as k — > 00, 
reflecting the relevant nature of the parameters u\^. 

The difference between the two channels is encoded only in different boundary conditions 
imposed on the solutions to this equation. In the case of the t channel, there is a bound 
state and we have chosen to expand in powers of the energy, po, relative to this state. The 
dimer propagator in this channel should therefore have a pole at po = in the physical limit, 
which leads to the condition 

u 1>t (k = 0) = 0. (17) 
The solution to Eq. ( TT6l) that satisfies this can be written in the form 



,,x 9 M g 2 1 
Ui,t{k) = — Tr77TT3 / d 1 



4vra f 2 (2 



7T 



1 1 



(18) 



g 2 M ( 4 2 k 2 1 . , 

^^ + ^73 7T + -COt" 1 (K)-l , (19) 



4ira t \3n 3tt (k 2 + 1) n 

where we have used Eq. (iTTj) for the deuteron energy and we have introduced the dimen- 
sionless variable k = kdt. 

The fact that SU(4) is not an exact symmetry means that the two scattering lengths are 
not equal. In particular, the s channel has no bound state and hence a negative scattering 
length, a s . If this were the only channel we were interested in, we could just expand energies 
around zero and the appropriate boundary condition would be 

u hs (k = = -\ , 20 

Ana, 



which is obtained by replacing a t by a s and setting Ed to zero in Eq. ( fl9l . However we 
cannot use this here since we need to treat both channels simultaneously and, as discussed 
above, we haven taken the deuteron energy as our expansion point. (A similar issue would 
arise even if the s channel also had a bound state but this was shallower than the deuteron.) 
The appropriate boundary condition can be obtained by replacing a t by a s in the first part 
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of the right-hand side of Eq. (fT9|) , while keeping Ed in the second term. Introducing the 
dimensionless parameter a (< 1) by 



l/a s = a/a t , 



we can write the solution in the form 



u 1)S (k) = ui )t (k) 



9 2 M 

4-KClt 



a 



(21) 



(22) 



This indicates that this approach is at least correct to first order in the strength of SU(4)- 
breaking, 1 — a. 

The wave function renormalisation factors both satisfy the same equation and we 
can impose the same boundary condition that they vanish as k — > oo on them. They then 
have the same form in both channels, 

11/"., 1 



z <f>A k ) 



4(2tt) 



d 6 q 



& + R N (q) - E d /2 

a t g 2 M 2 ( 2k (5k 2 + 3) 2 
- 2 -LI I 2 



S7T 



37T (K 



cot 1 (k) 



(23) 



Litim's cut-off 



20| respects Galilean invariance to the order we work here (see Ref. 21] 



for more details) and we find that the mass renormalisation factors are the same as those 
multiplying the energy, 

Z m ,i — Z<f>,i- (24) 

To simplify the expressions for the loop integrals in channels with more than two particles, 
we introduce a compact notation for the inverse propagators, 



E N R(q) 



q- 



2M 



£ d /2 + R N (q), 



AM 



+ w M (g)/Z - E d + R D (q)/Z (f> 



(25) 
(26) 



where we have suppressed the implicit k dependence of the running quantities in these 
expressions. 



B. Three-body couplings 

The evolution equations for the three-body (nucleon-dimer) couplings decouple into a set 
of four for the (1/2, 1/2) channel and two separate equations for the (3/2, 1/2) and (1/2, 3/2) 
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channels. They have the forms 

^(1/2,1/2) _ ^(1/2,1/2)^(1/2,1/2)^^ ^(1/2,1/2)^ 2 ^ 



+9 >4L 



x (1/2,1/2) x (1/2,1/2) 



Hs 



1 x (1/2,1/2) r 

2 A « i2 ' s 



+-<7 4 (J 3)t +9J 3) .) 
lb 



(27) 



a fc Ai 1/2 - 1/2) = Ai t 1/2 ' 1/2) AL 1/2 ' 1/2) /i,t + \i /2 ' 1/2) ^J 2 ' 1/2) h, 



+4^ 



o X (1/2,1/2) x (1/2,1/2) 



J 2 



3A (l/2,l/ 2) _ A (l/2,l/2) 



1 2 



a x (3/2,1/2) 



16 

x (3/2,l/2)\ 2 j 2 x (3/2,1/2) j 1 4 r 

A tt j J l,t + # A « L 2,t + -^9 1 3,t, 



(28) 
(29) 



where we have displayed only half the equations; the others can obtained by appropriate 
interchanges of spin and isospin labels, in particular t •f-)- s. We have also defined a short 
hand for the fc-dependant loop integrals 



hi 



(2tt) 3 
1 

(2^)3 
1 

(2^)3 



d 3 q 



QuRd + Zsd k R 



N 



(E DRti (q) + Z4,E NR (q)) 



2 ■ 



^3 E NR {q)dkRp + {E DRi + 2Z c f,E NR ) d k RN 
q ~ E NR (qy (E DRii (q) + Z^E NR {q)) 2 ' 

^3 E NR (q)d k R D + (2E DRi + 3Z ( f,E NR ) d k R N 
q ~ E NR (qr(E DRl (q) + Z^E NR (q)) 2 ~ 



(30) 
(31) 
(32) 



which can be evaluated in closed form for our choice of cut-off functions, Eqs. ( 12511261) . 



C. Dimer-dimer couplings 

The evolution equations for the four-body (dimer-dimer) couplings u 2 or u 2 also decouple, 
into a set of four for the (0, 0) channel and three separate equations for the (2, 0), (0, 2) and 
(1,1) channels. They have the forms (where again we have displayed only half the equations), 

d k u 2 , t = \u\ t K u - 2g 2 \f 2 ^K 2 - ^K 3 , (33) 



(34) 



d k u 2M = \u\ ts K Us - f\ U £ Af 2 ' 1/2) + £ A-j /2 ' 1/2)N j K 2 - ^K 3 , 

\ i=t,s ij=t,s / 

d k u 2tt = \u\ tt K u + X -u 2ts u 2st K^ s - \2g 2 \^ l,2) K 2 + ^K 3 , (35) 
d k u 2ts = ^u 2u u 2st K ht + -£u 2ss u 2ts K 1)S - 12g 2 \ { s t /2,1/2) K 2 - ^g 4 K 3 . (36) 
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Note that each equation contains only a limited subset of the thee-body couplings, as a 
result of the constraints of angular momentum and isospin. The loop integrals used here are 



K ^ = I d ^ 7 V > ( 37 ) 

(2tt) 3 J Z^EDR^q) 2 



Ki ts = j\-J d\ ^ Sj (38) 



f ,3 dkRjsr 



f 3 dkRN 



D. SU(4)-breaking form 

Since the two S'-wave scattering lengths for nucleons are very similar, SU(4) is a good 
approximate symmetry in nuclear physics at low energies. It is therefore sensible to re- 
express our couplings in terms of ones which multiply SU(4)-symmetric combinations of 
fields, and ones which are generated by SU(4)-breaking in the two-body sector. This is also 
need for analysing the behaviour of the couplings in the large-/c regime where we impose 
our initial conditions, as discussed above. Using the same notation as in Appendix A and 
denoting the SU(4)-breaking couplings with a "<5", we define, for example, 



A = \ (\^ l/2) + A£Ai/ 2 ) - 2X^) , (41) 

A' = ~ (AL 1/2 ' 1/2) + X^ + 2XT' m ) , (42) 
SX = X^-Xg^\ (43) 



"a = 2 ^ 2ss + Um + 2U2ts ^ ' ( 44 ) 

u 2 = 2 (^ss + u 2t t - 2u 2ts ) , (45) 

Su 2 = - {u 2tt - u 2ss ) , (46) 

for the (1/2, 1/2) and (0,0) channels. In a similar way we break the loop integrals up into, 
for example, 

h = \ (hs + ht) > (47) 

**i = ht - hs- (48) 
12 



The equations we obtain are generalisations of the SU(4) symmetric case discussed in 
Appendix A. For the spatially symmetric channels they become 

1 



d k X = hX 2 + g 2 I 2 X + -g 4 h 



d k SX 



d k u 2 



d k 5u 2 



~h6X 2 + 5X 



AT 

3 -5hX + 3 -g 2 SI 2 



d k X' = I 1 X ,2 -2g%X' + g% 



+- A h5X 2 + 5X 



3 -5hX> - \g 2 5I 2 



^Sh5X 2 + h5X(X + X') - ^g 2 I 2 5X 
+ -SI 1 XX' + g 2 -5I 2 {X' - 2X) - -g*5l 



6 

d k u 2 = \k x {u' 2 f - 2g 2 K 2 X - -^K, 
+ {&u 2 ) 2 + 5Ki5u 2 u' 2 , 



6' 



X -K x u\ - 2g 2 K 2 X' + ^g 4 K 3 



(5u 2 ) + 5Ki5u 2 u 2l 



-SKi (5u 2 ) 2 - 3g 2 K 2 5X 
+]-I<i5u 2 (u' 2 + u 2 ) + ]-5Kiu' 2 u 2 . 



(49) 



(50) 



(51) 



(52) 



(53) 



(54) 



Note that all SU(4)-breaking terms on the right-hand-sides are at least quadratic in sym- 
metry breaking quantities. We would thus expect that SU(4) is a rather good approximate 
symmetry in the three- and four-body systems. 



IV. RESULTS 

A. SU(4)-symmetric limit 

We first look at the case of exact supermultiplet symmetry, a s = a t . The evolution 
equations in this limit are described in detail in Appendix A. As explained there, we can 
reduce the problem to a limited set of parameters. These include two dimensionless three- 
body coupling constants, one of which (the coupling A' in the (1/2,1/2) channel) exhibits 



the Efimov effect 



12|. 



This effect is a remarkable feature of any three-body system with a attractive short-range 
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FIG. 1. An example of the full evolution of A' (solid blue line) compared to the asymptotic evolution 
(red dashed line). The full evolution was started by integrating downwards from t = 20, with the 
same initial condition as chosen for the asymptotic solution, [colour online] 

interaction. In the unitary limit (a —> oo) such a system possesses an infinite number of 
three-body bound states with a geometric spectrum. Even away from the unitary limit, 



these systems disp 
body observables 



ay universal features, such as relations between various three- and four- 
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231 ] . In the framework of a renormalisation gro up, the evolution of 



the corresponding three-body force shows a limit-cycle behaviour [19[, as consequence of 
the periodic appearance of new Efimov states as the cut-off is lowered. Away from the 
unitary limit, the finite inverse scattering length provides an infrared cut-off on the Efimov 
behaviour. This leaves a unique shallowest bound state which, in the nuclear context, we 
interpret as the triton. As a result, the periodic behaviour of the three-body coupling stops 
when the cut-off scale decreases to a value comparable with 1/a where it becomes almost 
independent of the running scale. More details of the running of A' in the unitary limit can 
be found in Appendix A. 

We first compare the asymptotic results for A', or rather a rescaled version of this coupling, 
A', defined as in Eq. (TT4"]) . to the numerical solution of the full equation. This coupling 
diverges periodically in t = ln(/«) = ln(ka t ), reflecting the appearance of the geometrically 
spaced states of the Efimov effect. In Fig. [I] we plot the arctangent of this coupling (after 
a suitable rescaling), which remains finite. The linear growth seen for large t is the signal 
of Efimov behaviour. The plot shows that the full solution follows the asymptotic form 
very closely down to t ~ 0. At this point, the scale 1/a* becomes important and acts as a 
low-energy cut-off on the tower of Efimov states. 

If we add the appropriate ie terms to impose causal boundary conditions on the prop- 
agators, we find that A' has an infinitesimal negative imaginary part. In our numerical 
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treatment we take advantage of this, by starting the evolution of A' at large k from the 
asymptotic solution, 

K 2 M(t) = ^31 - 5^535 tan ^ (5 - 5i - >/535*) ) , (55) 

with small a imaginary part 8. As we evolve downward in t, this allows our solution to 
bypass the singularities on the correct side and hence we are able to integrate the equations 
for the four-body couplings numerically. In this context we note that it is important not 
to take the finite imaginary part to be too small relative to the numerical precision used in 
the integration, as otherwise errors are introduced by the integration regions close to the 
singularities. Our results in the scaling region differ from those of Schmidt and Moroz {]], 
who took too small values for the imaginary part. 

The three-body parameter in spatially symmetric channel, A', couples to the evolution 
of the dimer-dimer parameter in the corresponding four-body channel, u^,. The rescaled 
version of this, which we call U2, is defined as in Eq. (I15j) . However in the region of large 
K ^> 1, it is more convenient to work with 

V dd (k) = k 3 U 2 . (56) 

The evolution of this in the scaling regime is shown in Figs. [2] and [31 There we can see 
that the solution, apart from a very short-lived transient, is completely determined by the 
driving term deriving from the coupling to A'. The coupling Vdd shows periodic singularities, 
reflecting the Efimov physics in the three-body channel. The initial transient effects vanish 
with a part of one Efimov cycle, reflecting the fact that the exact boundary condition on 
Vdd is an "irrelevant" parameter. More details of one of the singularities of Vdd are shown 
in Fig. [3j As we decrease the imaginary part of A', i.e., the value of 5 in Eq. ( l55i) . we see 
that the solution behaves much like a logarithmic singularity, with a sharp peak in its real 
part, and a finite jump in its imaginary part. 

The rescaled coupling Vdd goes to zero like k 3 as k — > 0. However we are ultimately 
interested in the values of the unsealed couplings in this limit. This suggests that the 
most appropriate technique for solving the differential equations numerically, is to use the 
equations for the rescaled quantities, Vdd and n 2 A' for t = In k > to ~ 0, and the equations 
for U2 and A' for t < t . Doing this we find that both couplings tend to finite values in the 
physical limit. The physical value of U2 is in general complex. This indicates that we are 
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FIG. 2. The evolution of Vdd as a function of t in the scaling regime, starting with the initial 
condition Vdd — at t = 100, for three different choices for the phase of A' on the asymptotic 
limit cycle (yellow dotted, red dashed, solid blue). The imaginary part of A' at the starting point 
is the same in the three cases shown, [colour online] 
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t t 

FIG. 3. A single singularity of Vdd in the scaling limit, for three decreasing values 5 = 10" 2 , 10" 4 , 
10 -6 (yellow dotted, red dashed, solid blue) of the imaginary part in the initial condition on A', 
Eq. (|55p . [colour online] 

indeed looking at an inelastic channel, as a result of the more-deeply bound states in the 
three-body channel. 

Finally, we express our results in terms of scattering lengths, as discussed in Appendix B. 
We take the value at = 4.32 fm for the two-body system in order to reproduce the deuteron 
binding energy. This differs slightly from the experimental value because our current cal- 
culations do not include finite range corrections (i.e., momentum-dependent interactions). 
Using the relations 

a DD /a t = 327rL7 2 (0), (57) 
a DN /a t = ^A'(O), (58) 

we get the results summarised in Fig. HJ These show the relationship between four-body 

1(3 
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FIG. 4. The relations between the scattering lengths in the SU(4) limit for the channels exhibiting 
the Efimov effect. All the points on these curves are obtained by choosing different points on the 
asymptotic limit cycle for A' as initial values, [colour online] 



scattering and the three-body parameter associated with the limit-cycle of A', which is fixed 
by the nucleon-dimer scattering length. There we see that the real part of the dimer-dimer 
scattering length goes through zero very close to the zero of the of the nucleon-dimer length. 
We also see that there is a maximum value for the real part of the DD scattering length. 

In the other pair of channels, which have spatial wave functions with mixed symmetry, 
all three- and four-body parameters are irrelevant. All the physical quantities can therefore 
be related to the single physical scale, at in our action. In particular, we find 



o'ddM = 327rt/ 2 (0) = 1.34, (59) 

e /215 

4 °\ if 29 
a' DN /a t = -A = -5L _ = -0.43. (60) 



The result that the dimer-nucleon scattering length in the mixed-symmetry channel has the 
opposite sign to the the nucleon-nucleon scattering length is a stable one. This is due to the 
fact that we have a quadratic equation for the coefficient k~2 in the asymptotic behaviour 
of A, with a discriminant b 2 — 4ac where a, b and c are all positive-and we thus have two 
negative solutions. There are thus two asymptotic solutions, both of which are negative. 
(Details can be found in the discussion surrounding Eqs. ( 1A24I) and ( 1A25I) .) The sign of 



A is preserved in the evolution out of the asymptotic regime, and hence, even though the 
numerical details can depend on the choice of cut-off function, as long as we insist on a 
consistent fc-scaling of the cutoff function the opposite sign remains. 
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B. Broken SU(4) 



In the real world, the singlet and triplet nucleon-nucleon scattering lengths are different 
and there is no exact SU(4) symmetry. The evolution equations do not decouple and so 
we have deal with 2x2 matrices of coupling constants in both the (1/2, 1/2) three-nucleon 
channels and the (0, 0) four-nucleon ones. In each case we can identify two scattering 
"eigenchannels" and, from the zero-energy T matrix, determine a scattering length in each 
of these channels. We show the general behaviour of these scattering lengths here, taking 
value of a = at/a s = —1/4 for the SU(4) breaking parameter. This is close to the realistic 
case, giving a s = —17.2 fm. 

The results depend on the starting on the limit cycle of the coupling constant that displays 
the Efimov effect in the scaling limit. The dependencies of the scattering lengths on this 
three-body parameter are presented in Fig. [5] Both lengths show avoided crossings; the 
narrowness of these is an indication that SU(4)-breaking effects are relatively weak. Further 
evidence of the smallness of this breaking is provided by the ratios between the components 
of the eigenchannel solutions. As can been seen from Fig. [6], the mixings are small (< 20%), 
except in narrow windows around the crossing points. 

In Fig. [7] we show the relations between the scattering lengths in the channels dominated 
by the Efimov effect. The scattering length in these channels vary rapidly with the three- 
body parameter. This means that at the avoided crossing, we switch from one branch 
of the solution to the other. As a result there is a very small discontinuity in the plots, 
close to the points where a on an d Refers] vanish. In the other channels we have an 
approximately constant scattering lengths away from the avoided crossing. These have 
values of a£,N — —0.26 ± 0.03 fm and Re[ao£>] = 5.55 ± 0.11 fm. 

C. Realistic scattering lengths 

We now turn to the implications of these results for realistic few-nucleon systems. Above 
we have already shown some results obtained with values for the nucleon-nucleon scattering 
lengths that are close to the observed ones. The small differences are not significant given 
the level of truncation of our REA, which means that effective-range corrections are not 
included. 



18 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0o 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



0o 0o 




0o 



FIG. 5. Evolution of the eigenvalues of the T-matrix (parametrised as scattering lengths) in the 
two B = 3 (1/2, 1/2) channels, and similar for the two B = 4 (0, 0) channels, as a function of the 
parameter 4>q specifying the initial condition on the limit cycle, [colour online] 



As already discussed, we also need one piece of three-body data to fix the starting point of 
the evolution on the Efimov cycle that is present for large cut-off scales. We choose to do this 



by using the experimental value of the spin-doublet nd scattering length 2 a^ d = 0.68 fm 24 1 
to determine the initial value of A^ 2 '^ 2 " 1 (since nd is not one of the scattering eigenchannels) . 
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FIG. 6. Ratio between components of the couplings in the eigenchannels for (left): A^ 1 / 2,1 / 2 ) and 
(right): u-i- [colour online] 



We have to be careful here, since, as stated earlier, we calculate the tree-body couplings off- 
shell, at energy Ed/2 below threshold. 

The situation is different for the nd scattering length 4 a n d in the spin-quartet channel. 
All the three-body parameters that contribute here are irrelevant and so the physical value 
of 4 a n d extracted at k = does not depend on their initial values. Our FRG calculations give 



■1.02 fm. This differs significantly from the experimental value of 4 a n d = 6.35 fm 



24j . Other theoretical calculations seem to be able to re pro duce this result, either by 31 1 



or by solving the Skornyakov-Ter-Martirosyan equation [32]. As we have said before, the 
negative value of this scattering length appears to be a stable result of our calculations, and 
is not very sensitive to the fine details. One possible reason for this is that the value of A a n d 
reported here is defined not at threshold but at the expansion point we have used for the 
nucleon energy in our flow equations, Ed/2~ — 1 MeV. Energy dependence is known to 
be important for three-body observables 27| and so that the extrapolation to the physical 
threshold could have a significant effect, especially since Ed is a small scale, we could expect 
significant energy dependence. The work of Bedaque et al [32| shows energy dependence, 
but when they include only the scattering length, this would probably not be sufficient to 
explain what happens here. In order to explore this in more detail, we will need to extend 
our REA to include energy/momentum-dependent couplings. 

In the case of four-nucleon systems, we consider two interacting deuterons and extract 
the dd scattering lengths in the spin-singlet and quintet channels. Besides being of interest 
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FIG. 7. The relations between the scattering lengths in the Efimov-dominated channels for a = 
— 1/4. [colour online] 

on its own, the low-energy dd interaction may have astrophysical applications, since under 
conditions expected to exist inside brown-dwarf stars a many-deuteron system may behave 



as a superfluid |28| . A framework like the FRG that can describe both the dd interaction in 
vacuum and be extended to dense matter would be very useful in this context. 

In the case of two deuterons, our expansion point corresponds to the physical threshold 
and so issues of extrapolation in energy do not arise. The singlet channel can couple to 
the n+ 3 He (or p+ 3 H) channel which has a lower threshold. Hence, as we have already seen 
above, the scattering can be inelastic, and its scattering length is complex. In contrast, 
the coupling of the quintet channel to the rearrangement ones is much smaller as non-zero 
orbital angular momentum is required 29[. In our S'-wave treatment, this channel is closed. 



Our result for v>2 corresponds to a singlet dd scattering length with Repa^] = 4.44 fm and 
Imparl = 0.17 fm. The real part is consistent with the value of Re[ 1 a^] = 4.9 fm obtained 



in Ref. 
is larger. 



30] by solving the Faddeev-Yakubovsky equation, even though the imaginary part 
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In the spin-quintet dd channel we get Re [ 5 add] = 2.55 fm. This also agrees with the 
value Repa^] = 3.2 fm obtained by Rupak in the framework of an e expansion 3lj. This is 
perhaps unsurprising since that work uses an effective field theory with a Lagrangian that 
like our REA omits effective-range terms. What is more puzzling is that Rupak finds a very 
different value from ours for 4 a n d = 4.78 fm, the scattering length in the three-body channel 
that feeds into the (2, 0) four-body one. This may be yet another suggestion that energy 
dependence of these couplings, allowing us to extrapolate better to the on-shell value, needs 
to be considered. 

Our result for the deuteron-deuteron scattering length in the quintet channel agree qual- 
itatively with the exact quantum mechanical analysis in Ref. 29[, given our incolplete treat- 
ment of the 3+1-particle rearrangement channels. The authors of that work find that this 
scattering length is very sensitive to these channels, and excluding them leads to a substan- 
tial reduction of Repaid], from 7.5 fm to —0.1 fm. 

which leads to a substantial reduction when the 3+1 break-up channels are excluded, 
although the value quoted in that paper is actually negative, Re [ 5 add] = —0.1 fm. 

In this first application of the FRG method to three- and four-nucleon systems, we are able 
to reproduce some of the effects of inelasticities in deuteron-deuteron scattering. However at 
this level we do not describe the nucleon-trimer threshold that is needed to get the correct 
energy dependence in this channel. To do this we would need to extend our approach by 
adding an auxiliary trimer field, along the lines suggested by Schmidt and Moroz This 
is certainly feasible but it would require adding a number of additional interaction terms to 
our REA, leading to a much larger set of coupled evolution equations. 



V. CONCLUSIONS 

In summary, we have applied the FRG method to three- and four-nucleon systems, both 
in the limit of exact Wigner SU(4) symmetry, and with realistic symmetry breaking. From 
the couplings in the physical limit, we have calculated the nucleon-deuteron and deuteron- 
deuteron scattering lengths in various spin-isospin channels. 

We find that the evolution of one three-body coupling shows oscillatory, limit-cycle be- 
haviour which is a manifestation of the Efimov effect in the corresponding channel. We 
therefore need to use one piece of data to fix one three-body parameter; all other three- and 
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four-nucleon observables in the spatially symmetric are then predicted in terms of this and 
the two-body scattering lengths. In contrast, the observables in the channels with mixed 
spatial symmetry are independent of the initial scale, provided it is chosen large enough. 
Their values are therefore determined by the two-body input alone. 

The Efimov effect appears in the spin-doublet nucleon-deuteron and the singlet deuteron- 
deuteron channels. Using the doublet scattering length to fix the three-body parameter, we 
get a value for the singlet deuteron-deuteron scattering length that is very close to one 
obtained in the exact quantum mechanical calculations. In contrast, our value for the spin- 
quintet scattering length shows no significant dependence on the three-body parameter, but 
differs from the results of other calculations, potentially as a result of ignoring off-shell effects 
in this work. 

There is a variety of ways in which the present study could be improved. One important 
one is the introduction of a trimer field in order to describe better the nucleon-trimer chan- 
nel, which is responsible for the inelasticity in deuteron-deuteron scattering. A second is to 
extend the ansatz for the running action to include energy- or momentum-dependent cou- 
plings. This would allow both the calculation of scattering observables away from thresholds, 
and the inclusion of effective-range corrections. Investigations of both of these are under 
way. 
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Appendix A: SU(4) limit 

In the limit of exact SU(4) symmetry, a t = a s , a = 1, the problem simplifies considerably. 
None of the loop integrals now depends on the type of dimer considered. Looking at the 



23 



A ) 


\ (V2,l/2) 
- A tt — 


\ (1/2,1/2) 


A') 


_ x (1/2,1/2) _ 
— A ts — 


X (1/2,1/2) 


A" 


_ , (3/2,1/2) _ 
— A tt — 


A (1/2,3/2) , 


u 2 


= U 2U + ^2fe, 




u 2 


= U 2t t - U 2 ts, 





evolution equations it becomes quickly obvious that it is useful to introduce new constants 

(Al) 

1 " " (A2) 

(A3) 
(A4) 
(A5) 

which satisfy much simpler equations than the "channel couplings". The ND coupling 
constants satisfy 

d k X = I 1 X 2 + g 2 I 2 X + ^g%, (A6) 
d k X' = I 1 \' 2 -2g 2 I 2 \' + g%, (A7) 
d k X" = h\" 2 + g 2 I 2 X" + \g%. (A8) 

Since the first and last equation are identical, we thus realise that A = A". The potential 
then simplifies to 

V 3 ,BF = A' [t t (p30,P 3 )^ t (PlO,Pl) - S t (P30,P 3 )r(PlO,Pl)] 

[*P(P20,P 2 )t(p40,P4) ~ ^(P20,P 2 ) S (P40,P 4 )] (1/2 ' 1/2) 

(1 /2 1 /2) 

[^(P20,P 2 )^(P40,P 4 ) + ^(P20,P 2 ) s (P40,P 4 )] (1/2 ' 1/2) 

+A [s ] (p3o,P 3 )^(pio, Pl )f /2 ' 1/2) ■ mP2o,P 2 )s(p i0 ,P,)f M2) ■ (A9) 
The u 2 equations also simplify 

d k u 2 = X -u\K x - 2g 2 \"K 2 - ^g 4 K 3 , (A10) 

d k u 2M = \u\ M K x - 2g 2 XK 2 - ^g 4 K 3 , (All) 

d k u' 2 = U' 2 ^ - 2g 2 XK 2 - ^g 4 K 3 , (A12) 

d k u 2 = X -u\K x - 2g 2 X'K 2 + ^g*K 3 . (A13) 
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Thus we also conclude that U2,t s = u 2 = u' 2 . The potential terms in this limit can thus be 
written as 



a simple sum of the non-local generalisations of the square of the linear and the quadratic 
Casimir invariant of U(A). 

1. Unitary limit 

It is instructive to look at what happens for large k = ka^. This corresponds to the 
unitary limit, where the scattering length diverges. The behaviour for large k is universal, 
and thus provides boundary conditions for the numerical solution of the full problem. 

a. Three-body coupling 

In the three-body case we are left with a set of dimensionless integrals (the tilde denotes 
we have scaled the expressions give before by appropriate powers of M and a t ) 




( ^(PMjPa) ■ *(P40,P 4 ) + S G°20,P 2 ) ■ S G°40,P 4 )J 

+ Y^2 - U 2) [* + (PlO,Pl) ■ t\p 2 Q,P 2 ) - S f blO,Pl) • S t (P20,P 2 )] 
[*(P30,P 3 ) • *(P40,P 4 ) - S (P30,P 3 ) ■ S(P40,P 4 )] , 



(A14) 



16/t 2 (k 2 + 1) (45 (k 2 + l) 3 cot _1 (K) + k (67k 4 + 120k 2 + 45) 



i 



5 (3 (3k 2 + 10) (k 2 + l) 2 cot -1 («) + k (31k 4 + 59k 2 + 30)) 2 



(A15) 




32k 4 (15 (6k 2 + 13) (k 2 + l) 2 cot _1 («) + k (222k 4 + 415k 2 + 195) 



(A16) 



(A17) 



5 (k 2 + 1) (3 (3k 2 + 10) (k 2 + l) 2 cot _1 (K) + k (31k 4 + 59k 2 + 30)) 2 



The differential equations now read 




(A18) 



(A19) 
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FIG. 8. Four examples of the evolution of k 2 A. These were solved starting the integration at 
t = In k = 20. The three dashed lines represent the solutions for three random choices of the initial 
value. The solid (green) line starts using Eq. (|A25p as initial condition, [colour online] 



The 7's have simple asymptotic expansions for large k, 

h = k \ 
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The asymptotic solution to the equation for A flA18j) can be obtained using these results as 

A(k). 
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2 rmr 
CAtsV 5 



(A23) 



It is very tempting to assume c ^ 0, and get the asymptotic solution for k — > 00 



(A24) 



Numerical investigations, Fig. [HI of the solution to the "full" equation, where we have not 
made the asymptotic expansion, starting the integration from t = ln(ft) = 20, shows that 
when we have chosen an arbitrary initial condition at finite k, the solution typically has a 
quick decay to the fixed point value for c = 0, for almost all initial conditions. Only when we 
start very close to the asymptotic value (1A24I) do we find an apparently divergent numerical 
solution. Thus it is actually most efficient to use the asymptotic solution with the choice 
c = 0, 
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(A25) 
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as initial condition for large k. 

The differential equation for A' does not have such a trivial fixed point; the asymptotic 
solution is found to be quasi-periodic 

k 2 A'(k) = ^ ^31 + Sv^tan ^v / 5351n(fi;) + <fo^ . (A26) 

This limit cycle-behaviour is the basic manifestation of the Efimov effect in the approach 
employed here. 



b. Four-body couplings 

For the fur body-couplings we can also define scaled integrals, 

g 192tt 2 k 4 (k 2 + if 

1 ~ 5 (3 (k 2 + l) 2 cot _1 («) + k (5k 2 + 3)) 

f 15 (k 2 + l) 3 cot _1 (K) + k (17k 4 + 40k 2 + 15)' 

X 2 » ( A27 ) 

(3 (k 2 + l) 2 (k 2 + 8) cot _1 (K) + k (21k 4 + 43k 2 + 24)) 
4k 4 

K 2 = (A28) 

3tt 2 (k 2 + 1) 3 V ; 



Ks = ^77-7-^4- (A29) 



8k 

3tt 2 (k2 + 1. 
We can rewrite the differential equations as 



d K U 2 = - 2A'A 2 + ^K 3 , (A30) 

d K U 2 = \\Jlk x - 2AK 2 - ^K 3: (A31) 



and the asymptotic expansions 



47T 2 K 2 

K X = (A32) 
K 2 = -4^, (A33) 



37T 2 K 2 

8 

3tt 2 k 4 



As = (A34) 



It pays off to use 



V dd (k) = k a U 2 (A35) 



This satisfies 



Kd K V DD = %-V dd {k) 2 + 3V dd (k) + \- -^k 2 A'(k). (A36) 
15 n z o% z 
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Without the additional A term we find for large k that 



V DD _« + v^545 m _ 2 13551 

47T 

By numerical solution we find that the periodic oscillations in A dominate over this be- 
haviour. 



2. Scattering lengths 



a. Dimer-dimer scattering 



Our effective action contains the term 

V DD =-U 2 J S {4) {p 1 +p 3 -p2-p4)[(pUPw,P 1 )<Pa{P20,P2)] $(P30, P 3 )^(P40, P 4 ) 



+ 



/ <^ (4) (Pl +P2 -P3 -P4) [* f (PlO,Pl) • t\p 20 ,P 2 ) ~ S f (PlO,Pl) ' S j (p 20 ,p 2 



X [*(P30,P 3 ) ■ *(P40,P 4 ) - S (P30,P 3 ) ■ S (P40,P 4 )] 



(A38) 



This is in terms of classical fields; the corresponding operator in the "pole approximation" 
when the cut-off reaches zero is 



V D d =t; u 2- n 2 



t\E B , 0) • t\E B , 0) - s\E B , 0) ■ s\E B , 0) 
[t(E B , 0) • t(E Bl 0) - s(E B , 0) • s(E B , 



— u 2 
12 



(A39) 



--u' 2 n (n - 1) 
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"2 



t f (£ B , 0) • £ f (£ B , 0) - ^(Eb, 0) ■ s\E B , 0) 



[t(£7 fl , 0) • t(E B , 0) - s(E B , 0) ■ 



(A40) 
(A41) 



Note the important effect of normal ordering; this is simply a re-statement of the fact that u' 2 
does not contribute in the single dimer channels. We now calculate the scattering between 
the normalised scalar ([2] denotes the symmetric irrep of SU(4) [3^|) states 



l[2](0,0)> 



^6x2 



and we get 



t\E B , 0) ■ t\E B , 0) - s\E B , 0) ■ s\E B , 0) |0) 



([2](0,0)|V D b|[2](0,0))=u 2 . 



(A42) 



(A43) 
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We thus find that 



M B _ 2 M N _ s 

aDD = 17 U2/ * = ~27 U2/ 4 

= a t 32ir(U 2 ) 



(A44) 



We also use the fact that Z^O) = a t g 2 M 2 /8ir. 

In the same vein we find in the other channel ([1,1] denotes the antisymmetric irrep of 
SU(4)) 

t\E B , 0) • t\E B , 0) + s\E B , 0) • s\E B , 0)1 |0) , 



|[M](o,o)> 



(A45) 



that 



M 



AS _ Mn /r7 2 

l DD 



a t 327r(U 2 ) 



(A46) 



In the mean-field limit, we find that lim^o {7 2 = 1/167T, and we get the standard result 
a r> S D = ^ a t- 



b. Dimer-nucleon scattering 

There are two possible dimer-nucleon state we can construct, a symmetric and antisym- 
metric one. 

We first assume that the scattering is evaluated in the channel of the symmetric dimer- 
nucleon states 



\DN) = 



iHE B ,0)^(E B /2,0)-s<(E B ,0),pHE B /2,0) |0), (A47) 



1/2,1/2 



J ms,m s 



we get 



(DF\i\DF) = A. 



(A48) 



We thus conclude 



M red 2i4/3. /7 

a DN — — — — A/Zj^ 



Air 

4 A 
at- A. 
3 



(A49) 
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In the antisymmetric case, defined by 



1 x 1 1/2,1/2 

\DN') = — &(Eb, 0)4>\E b /2, 0) + s\Eb, 0)ip\E B /2, 0) |0> , (A50) 



V2 

we get 

M pedx ,._ 2Mjv/3 v 



m s ,m s 



(I 



-X'/Z^ - — X'/Zj 



v ""4vr 4tt 
= a t ^M. (A51) 



Appendix B: Link between T and V 

Our main results are for the threshold behaviour of various quantities are, in the main, 
expressed in terms of scattering lengths. The relation between our effective action and a 
scattering length is based on the definition 

HmT(fc) = —at. (Bl) 

k^O M red 

For identical mass, Mi = M 2 = M, we get 

T(0) = ^a t . (B2) 

What we call a "potential" in the effective action, is really nothing but the classical form 
of the T matrix. For dimers we find a similar relation, but we have to interpret T as 
the expectation value of the quantised potential in the relevant channel. As shown in the 
appendix, this allows us to take into account symmetry aspects of the problem. Thus 

V = (V) = T DD (0) = (B3) 

which can be used to show that the dimer-dimer scattering length in the mean-field limit 
goes to 2a t . 
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